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1 Introduction 

This note introduces a PDE-like Toy-Model that embeds several aspects of Territory Working. The 
long term goal is to build a software tool that behaves like a Territory, and in particular that incor- 
porates its multi-scale-in-time-and-space nature, in order to make simulations, to explore scenarios, 
to foresee policy impacts, to help to make a decision when facing a change in the environment or 
in cultural behavior, etc.. The aim of this note is prove that the concept of building a model that 
couple Systemic Approach and PDE tools to achieve the evoked long term goal is possible. 
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Figure 1: Systemic Approach of Territory Working 



2 Brief Systemic Approach of Territory Working 

By Territory, we mean a Country or a Town or a set of Towns and Countries having a coherence 
and being administrated by a common local government. 
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The Woking of a Territory results from the interactions between several compartments. A 
representation of this Systemic Approach is given in Figure [T] Three compartments are considered. 
The first one ("People") concerns the population, the second one ("Energy") concerns the energy 
questions and the last one ("Economy") concerns the world of enterprises . Then, the compartments 
influence each other or themselves. This fact is symbolized by the arrows. For instance, the arrow 
which is the more at the top of the figure translates that when people are at work, they contribute 
to the production of wealth. In return, the arrow just below means that displacements are induced 
by enterprise locations and the arrow in the left translates that an energy consumption is induced 
by those displacements. The three arrows that point on the "Economy" compartment from itself 
symbolize that enterprises create wealth and that this created wealth is used by the enterprises 
for growing up and for efficiency improvement. The two last arrows, in the bottom, express that 
enterprises consume energy to create wealth. 
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Figure 2: Initial Population Density Vq (left) and Job-station Density £q (right) 



3 The Toy-Model 



In this Toy-Model, all quantities are dimensionless. The simulations are done using |Freefem- 



All the simulations are done over the same time interval which is run by variable t G [0, T) for a 
real number T > 0. Yet as the characteristic times of the various phenomena are not the same, 
they are given with various steps. 



The geographic model of the Toy-Model is a disk D, provided with coordinates x, with a boundary 
9D. The vector of norm 1, orthogonal to 9D and pointing outside D is denoted ji. 
On this disk, several densities, which depend on time, are defined: 



• V = V(t^x) is the Population Density. 

• 8 = S{t^x) is the Job-station Density. 
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Figure 3: Job-station Density Sq times Efficiency Indicator ig 



• W = W(t, x) is the Wealth-that-goes-to-people Density. 

For those densities, initial values are defined: Vo = Vo{x), which is drawn on the left of Figure 
[2I and which is the Population Density when everybody is at home at night; Sq = £o{x) which is 
drawn on the right of Figure [2I and, Wo = yVo{x) = 1. 

An Efficiency Indicator i = i(t, x) is also defined on the disk. It ranges in [0, 1] and mea- 
sures the efficiency of the enterprises at every point of the disk; it is initialized at Xq. The product 
of the initial Job-station Density £0 by the initial Efficiency Indicator to = io(^) is given in Figure [3J 

From initial Population Density Vo and Job-station Density S two attractors are built. The 
first one As = As{t^x) attracts the Population Density towards Job-stations. It is defined as 

/ Vo{x)dx 

As{t,x) = — ^ {Vo{x)^a£{t,x)) , (1) 

/ Vo{x) dx -\- a £{t^ x) dx 

where a is a constant larger than 1. The second one Aj> = Av{x) makes the population to go back 
home. It is defined as 

Av{x)=Vo{x). (2) 

With those two attractors, we can simulate the daily motion of the population. For this, we 
firstly adopt the convention that day length is 1 and we define two functions tm = im{t) and 
te = te(^) depending only on time and periodic of period 1. t^(t) > for instants corresponding 
to the morning displacements, when people go to work and otherwise; and, ie{t) > for instants 
corresponding to the evening displacements, when people go back home and otherwise. Secondly 
we write the equation meaning that in the morning the Population Density is attracted by Job- 
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Figure 4: Displacements from home to work and back (Animation available on YouTube: |http : 
7/youtu . be/LkPlVT- a4pg[) 



station Attractor and in the evening by Home Attractor: 
dV , 



dt 



(t, x) = im{t) {As{t, x) - V(t, x)) + ie(t) {Av{x) - V(t, x)) , Vx G D, Vt G (0, +00), (3) 



The evolution of V is given by the movie in Figure [4j In this movie, we can see the alternation of 
people displacements to their works (in the morning) and to their homes (in the evening). Because 
of Job-station distribution most of people converge everyday to work in a small region located near 
the center of the disk. 



This motion generates an energy consumption. It is modeled by a time density (j) = (j){t)^ which 
is constant per day and which value at any given day is the double integral over the disk of the 
distances from home locations to Job-station locations weighted by the locations where Population 
Density V increases and where it decreases. 
In other words, ^ is a constant over every interval [n, n + 1], n G N with worth 



b-= f J\X- yfVrncrix)V^ecM dxdy, 



(4) 



where T^i^^r ^^ ^^^ density of increasing population, where population increases, when comparing 
the Population Density in the morning (before people go to work) and in the middle of the day 
(when workers are at their job-station). It is defined by 



^rncr(^) = max(7^(n + ^ ' ^^ ~ ^^^' ^)' ^)- 



(5) 



In a similar way, "Pgecr i^ ^^^ density of decreasing population, where population decreases, when 
comparing the Population Density in the morning and in the middle of the day and is defined by 



PB,„(a;) = max(-P(n +^,x) + V{n, x), 0). 



(6) 
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Figure 5: Evolution of Job-station Density £. (Animation available on YouTube: http : //youtu . 
be7PorgK 4q5Ds0i) 



The fact that people is at work produces Wealth with a rate which depends on the product of 
the Population Density V by the Job-station Density £ times the Efficiency Indicator i. To model 
this, we introduce a wealth production rate density uo = uo{t^ x) which is constant per day and which 
value at any given day, represented by interval [n, n + 1], n G N, is 



cj^(x)=/3o7^(n+-,x)£:(n+-,x)i(n+-,x). 



(7) 



where Pq is a coefficient small in front of 1 , meaning that the time scale of Wealth variation is large 
when compared to a day. 

This Wealth is split into three parts. The first one is allocated to enterprises' growth. In other 
words, it is a source term of the ODE, set in every point of the disk, the Job-station Density 8 is 
solution to: 



d£_ 



{t,x) = f3iuj{t,x), yx G D,Vt G (0,+oo), 



(8) 



where Pi is a coefficient belonging to a triplet (/3i, 1^2 •, P^) of positive coefficients such that /3i +/32 + 

The second part (/?2w(t, x)) contributes to improve enterprises' efficiency. It contributes as a 
source term in the ODE, set in every point of the disk and which is such that its solution i always 
ranges in [0, 1]: 



Si 
di 



{t,x) = \{t,x){l - \{t,x)) {I320j{t,x) - I3i0j{t,x)) , Va;eD,Vt e (0,+oo). 



(9) 



The term Piuoit^x) in this equation models that when organization grows, if nothing is done, its 
efficiency decreases. 

The result of the influence of Wealth on the Job-station Density £ and on the product of the 
Job-station Density £ times the Efficiency Indicator i are given in the movies of Figures |5] and [6J 



We see in those movies that, around the smah region where most of people converges everyday to 
work, there is a strong increase of Job-station Density £ and of enterprises' Efficiency i. 
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Figure 6: Evolution of Job- station Density £ times Efficiency Indicator i. (Animation available on 
YouTube: |http : //youtu . be/B7g jmu7q8iul ) 



The third part of Wealth production will diffuse inside the population. It is then a source term 
in the diffusion equation W is solution to: 

— -(t,x) -y/\yV{t,x) =/33Cj(t,x) + A^(^(t), VxGD,Vt G (0,+oo), (10) 

—- = J^w, VxGaD,VtG (0,+(X)), (11) 

where v is the diffusion coefficient of Wealth within the population on the territory. In this PDE, 
the energy consumption time density (j) is also a source term but with an opposite action than that 
of the Wealth production. This is because (j) induces a consumption of Wealth. In this PDE there 
is also a boundary condition on the disk border, translating the fact that there is a little wealth 
that enters the Territory through its border. The evolution of Wealth-that-goes-to-people Density 
is given by the movie in Figure [7[ We can see the diffusion of the Wealth-that-goes-to-people from 
the small region where most of people converges everyday to work towards the other parts of the 
Territory. 

4 Interpretation 

Clearly, this Toy-Model does not pretend to give realistic results. Nevertheless it seems that the 
proposed approach brings a way to couple several aspects of Territory Working. 

For instance, it can handle non-linearities. This is illustrated by the two pictures of Figure [8J 
The one on the left is the Population Density at midday of the first day and the one on the right is 
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Figure 7: Evolution of Wealth Density W. (Animation available on YouTube: http : //youtu . be/ 
OSJNyznPxvwp 



the Population Density at midday of the last day. We can see that, in the small region close to the 
center of the disk, the density is higher in the right picture than in the left picture. This can be 
explained as follow: Since people is working their, wealth is produced, which induces Job-station 
number to increase. As a consequence, day after day, more and more people is coming to work in 
this small region. 



The scales (regarding time and space) have no realistic meaning. Nevertheless, we can see the 
capability of this kind of models to account for a wide variety of scales. For instance, regarding the 
space scales, in the Toy-model, there is the size of the Territory (the disk) and there are the char- 
acteristic sizes of variation of Population-at-home Density, of Job-station Density and Enterprise 
Efficiency. Beside this, this way of modeling allows us to connect several Territories through their 
common boundaries and then to consider a large number of connected Territories. 
Regarding the time scales, in the Toy-Model there are the characteristic time of Population motion, 
the characteristic time of wealth production and the one of wealth diffusion. 
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Figure 8: Population Density at midday of the first day (left) and of the last day (right) 



